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We point out that current calculations of inclusive two-particle correlations in p-A collisions 
based on the Color Glass Condensate approach exhibit a contribution from Hanbury-Brown-Twiss 
correlations. These HBT correlations are quite distinct from the standard ones, in that they are 
apparent for particles widely separated in rapidity. The transverse size of the emitter which is 
reflected in these correlations is the gluonic size of the proton. This raises an interesting possibility 
of measuring the proton size directly by the HBT effect of particle pairs produced in p-A collisions. 


I. INTRODUCTION 

The observation by the CMS [1] and ATLAS [2] collaborations of ridge correlations in high multiplicity p-p collisions 
at the Large Hadron Collider (LHC) has triggered a wave of theoretical attempts to understand the nature of such 
correlations. Similar correlations have been subsequently observed by all three large LHC experiments in p-Pb 
collisions [3], and much more detailed studies of the properties of these correlations are available today. 

There are indications that the origin of these correlations is due to collective (hydrodynamic?) behaviour of 
the system produced in the collision [4]. However, a good quantitative description of the data is also achieved [5] 
within the Color Glass Condensate (CGC) based “glasma graph” approach [6, 7], which ascribes the origin of the 
correlations entirely to the structure of the initial state. In reality it is likely that both mechanisms are contributing 
to the correlations, probably in different transverse momentum ranges. At any rate, it is interesting to have a deeper 
understanding of the initial state induced correlations. Note that other explanations within the CGC [8, 9] and in 
other frameworks [10] also exist, but they will not be touched upon in this paper. 

In a recent paper [11] we have shown that the bulk of the correlations in the glasma graph approach originates from 
the Bose enhancement of gluons in the incoming projectile and target wave functions. This quintessential quantum 
mechanical effect increases the probability to find gluons with the same transverse momentum in the wave function 
of the incoming projectile (and target). The scattering smears the momentum of the gluons to some degree, but 
significant correlations in the gluons produced in the final state are nevertheless observable in favourable kinematics. 
In this note we point out to an additional physical effect present in the glasma graph calculation: Hanbury-Brown- 
Twiss (HBT) correlations between gluons far separated in rapidity (for reviews on interferometry studies in heavy-ion 
collisions see [12]). This leads to a potentially observable effect in the final state mesons which may allow a direct 
measurement of gluonic “size” of the proton. The HBT signal is also interesting in that it correlates gluons with same 
and opposite transverse momenta (in analogy with the double ridge structure of [5]). 

The fact that the CGC approach contains HBT correlations is not new and has been recognised before [13, 14]. 
The aim of this note is to understand their unique features, in particular their long range in rapidity nature and the 
fact that they reflect the gluonic size of the proton, within the standard setup in HBT studies [15]. We also make 
contact with the glasma graph approach. 

This paper is organised as follows. In Sec. H we present a short review of HBT correlations. Here we essentially 
follow the excellent reviews [15] with some slight change of notation. In Sec. HI we discuss HBT correlations in a 
system of gluons emitted from a Lorentz contracted source, and point out that in this case the HBT correlations 
extend over a long range in rapidity. We also discuss the typical transverse structure of the emitter expected in the 
CGC approach, and note that in p-A collisions its transverse size is that of the proton E?, while the transverse area 
of the region from which the emission is coherent is the inverse saturation momentum of the target Thus, if 

^ 1 the emission is dominated by the HBT signal that measures the proton size. In Sec. IV we translate 
this qualitative discussion into the language of glasma graphs and show which glasma graphs correspond to the HBT 
signal discussed above. In Sec. V we offer some concluding comments. 
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II. BASICS OF HBT 


In this section we review the basics of HBT correlations, following closely the reviews [15]. The object of study is 
the normalized two-particle correlation function. The single and two-particle spectra are defined as 


„ , , dN ^ , 
^i(P) = ^ = 


^2(Pa,P&) = 


dN 


d?Pad^Pb 




( 2 . 1 ) 

( 2 . 2 ) 


where and d are the creation and annihilation operators of particles (usually charged mesons) with momenta 
Pi = (p°, Pi) and (• • •) stands for the averaging over the ensemble of events. 

The two-particle correlation function C'(pa,p&) is defined as 


C'(Pa,Pb) = 


{NY 


^2(Pa,P6) 


{N{N -l))Vl{pa)Vl{pb) 
where {N) is the average number of particles and {N{N — 1)) the average number of pairs, such that 


I d^pPYp) = {N), 

J d^Pad^PbP2{Pa,Pb) = {N{N - 1)) . 


(2.3) 

(2.4) 

(2.5) 


The averaging over the event ensemble can be represented in terms of a density matrix p, so that for an arbitrary 
observable O, one has 


(6)=tr(pd). (2.6) 

The particles are emitted from the interaction region, and stop interacting once they leave the “freeze-out surface”. 
Once outside the interaction region, the particles propagate freely until reaching the detector. The standard description 
of the “emitter” is in terms of a classical source J{x) which emits the pions, so that for a given source J{x) the state 
of the pion field (between the emitter and the detector) is given by 

|J) = e”/^exp J d^pJ{p)dl^ |0), (2.7) 

where J(p) is the on-shell Fourier transform of the classical source J{x), 

/ d^ X 

—=^= exp[i{Ept-p-x)]J{x), (2.8) 

v 

with Ep = 'Yp^ + and the normalization of the state n = f d?p | J(p)p. This coherent state is an eigenstate of 
the annihilation operator. 


apjj) = iJ{p)\J) . 


(2.9) 


In the case of the emission of the particle from a single coherent source J{x), the density operator of the ensemble 
is just a projection operator on the coherent state: 


p=|J)(J|. (2.10) 

One then obtains 

Pi(p) = (J|atap|J) = |J(p)p, (2.11) 

^2(Pa,Pb) = ('/|ap„ap,ap„apj>/) = [•/(Pa)!^! J(P6)r (2.12) 


The two-particle spectrum is given by the square of the one particle spectrum and thus there is no Bose-Einstein 
correlations from the emission off a single coherent source. 
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A fixed source J{x) corresponds to a single event. Averaging over ensemble of events corresponds to averaging over 
the ensemble of sources J. In reality different regions of the freeze-out surface are not coherent. They vary from 
event to event independently of each other. This situation is modelled by considering the emitter to be a chaotic 
superposition of classical sources, whose phases vary in a random uncorrelated way. As a model for such source we 
take 


N 

J{x) = - Xi). (2.13) 

i=l 

Here each individual source has been shifted from its original (four-)position by Xi and boosted by momentum pi. 
This of course corresponds to the individual sources to be spatially separated as well as moving independently of each 
other. Additionally, a random phase ipi is associated with each individual source Jq{x — Xi). The total emitter J{x) 
is the sum of all the individual sources. The on-shell Fourier transform of it is written 

N 

= (2.14) 

2=1 


where 


Mp-Pi) = 


drx 




Jo(a:). 


(2.15) 


The state of the produced particles now depends on the source positions, source momenta and phases. The ensemble 
of sources, corresponding to the distribution of events, can be defined in terms of a density operator p which specifies 
the distribution of the source parameters. Assuming that the number of sources is distributed with a probability Pat, 
the phases are distributed randomly between 0 and 27r, and the source positions Xi and momenta pi are distributed 
with a density n{x,p) with the following normalization: 


^ Pat = 1, NPn = Ns, d'^x d^pn{x,p) = 1, 


iV=0 iV=0 

the corresponding ensemble average of an arbitrary operator O is given by 


tr(/50) 


N^O 


Dn[J] {J[N; {x,p, (l)}]\d\J[N; {x,p, (j )}]), 


=0 


where the measure D]\[[J] is defined as 


N 


Dn[J] =Y\_d‘^x^d‘^Pin-ixi^Pi)^ ■ 

i=l 

Using Eq. (2.17), it is straightforward to calculate the single particle spectrum: 

oo „ „2Tr JV 

(atap)=^Pjv/ / Dm[J] Y 


N^O 


(2.16) 


(2.17) 


(2.18) 


(2.19) 


=0 


Due to the averaging over the phases only the terms n = n' give non-vanishing contribution. Hence, the single particle 
spectrum reads 


'Pi{p) = Ns / d'^x'(Tp'n{x',p’)\Jo{p-p')\ = {\J{p)\ ). 


( 2 . 20 ) 


Before calculating the two-particle spectrum let us consider 

.277 N 


oo 


Nl I Dn[J] Y 

N=0 " n,n' = l 


Ar=o 

= Nsj d^x'dyn{x',p')e^^P>’-P^^-^'j*ipa-p')Mpb-p') = {J*‘^iPa)jHpb)). 


( 2 . 21 ) 
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Using Eq. (2.21), we get for the two-particle spectrum 




a®pi>) “ ^ ' Pn 


N=0 


N 


dn[j] 

n,n' — l 


g j(0n-l-0m-0„/-0m') g*Po ■ -a;„/) gipb ■ (Xm-a;^/) 


X Joipa -Pn')JoiPb -Pm')Jo{Pb -Pm)JoiPa - Pn) ■ 


( 2 . 22 ) 


The phase averaging leaves three non-vanishing contributions: n = n',m = m'; n = m',m = n' and finally the 
contribution which corresponds to emission of particles from the same source Jqj namely n = n' = m = m!. The first 
two contributions factorize a la Wick into products of the single particle averages similar to Eq. (2.21). The third 
contribution cannot be written in this form, and we will call it “irreducible”. It is affected by the structure of the 
current Jo- In fact, the coherent state approximation to an individual emitter Eq. (2.7) is not always appropriate. 
Corrections to this classical structure of the source do not affect the factorizable contributions in Eq. (2.22), as 
those contributions depend only on the single particle spectrum of an individual emitter. They do however affect the 
irreducible contribution, since it probes possible correlations within a single coherent emitter Jq. However when the 
average number of coherent emitters Ns is large, the irreducible contribution is suppressed by the factor 1/Ns- We 
will come back to this point later in Sec. IV. 

All in all, the two-particle inclusive spectrum can be written as 


^®Pa®P6®P“®Pi>^ — 'y y / I IT ^ Xi d Pin{Xi,Pi) > ^ ^ S \ JoijPa Pn)\ \JoiPb Pm)\ 

N=0 J n/m ^ 

+gdP.-P.)-Un-Xm) - Pm)JS{Pb - Pn)MPb - PmVoiPa " P„)| + , (2.23) 


where the is defined as the n = n' = m = m' contribution in Eq. (2.22). Thus, the two-particle spectrum 

reads 

^2(Pa,P.) = ^^^^^~^^^ {(|i(Pa)P)(|J(p.)P) + |(J*(Pa)J(p.))p} + (2-24) 

where we have indicated explicitly the 1/Ns suppression factor in front of the irreducible contribution. Neglecting for 
the moment the irreducible contribution, the single and two-particle distributions can be written in terms of the so 
called ’’emission function” S{x,K) defined as 

and closely related to the Wigner distribution of the produced particles [16]. The correlators of the source functions 
that appear in both single and two-particle distributions are written in momentum space. Eourier transforming them, 
we write 


Setting X = \{xi + X 2 ), y = xi - X 2 , q = Pa - Pb and K = ^{pa +Pb), we get 

Hence, the two-particle correlation function can be written as 


(2.26) 


(2.27) 


C'(q,K) = l 


\Jd‘^xS{x,K) 


/ d'^xS [x, K -I- iq) J d'^xS {x, K — ig) 


(2.28) 


The second term in this expression expresses the HBT correlations, which arise due to the large number of incoherent 
sources that constitute the emitter. 

The qualitative features of this expression are easy to understand. As a simple example let us take the individual 
sources to be static and distributed inside some radius R. For simplicity we take the actual distribution to be Gaussian, 
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n{xi) oc exp{ —Also for simplicity we will assume that each source has a Gaussian profile with a radius a much 

2 

smaller than R, Jo{x) oc exp{ —^}. For the (square root of) the numerator of Eq. (2.28) we obtain 


J d^XiC ^ j Xi+l)^] 


= I d^xd^yd^XiC 2a^ e 


(2.29) 


As long as a <C i?, the second Gaussian acts as a delta function, not letting Xi to stray too far from x, and thus the 
integral over Xi gives 


d^xd^ye 4 ^ oc e ’2 , 




(2.30) 


In this simple model the HBT signal (the correlation between particles at different momenta) has a width ^ 1 /R in 
the momentum difference which does not depend on the total momentum of the pair K : 

G(q,K) = l + . (2.31) 

In reality the shape of the emitter is more complicated and has to be modelled in some way. Also parts of the 
freeze-out surface are moving, and the source only exists for a finite time. All this introduces dependences on K as 
well as on the distribution of pi as well. However, for our purposes the simple example above suffices since it illustrates 
the basic physics of the HBT correlations. 

So far this has been a review of the standard theory of HBT. Our goal, though, is to understand what kind of 
HBT correlations one expects in the CGG picture of the collision. There are several important aspects which set this 
picture apart from the one that we have just described. 

First, one is interested in the emission of gluons and correlations between emitted gluons. For the sake of the 
argument we are going to assume local parton-hadron duality and forget about hadronization corrections. 

Second, in the GGC picture the collision is boost-invariant. The colliding objects are strongly Lorentz contracted, 
and they overlap with each other only for an infinitesimally short time. The emitter exists only for this very short 
time, and this has a profound effect on the nature of the HBT correlations. 

Additionally, since we are discussing p-A scattering, we assume that the saturation momentum in the nucleus is 
much larger than that of the proton. Recall that the inverse saturation momentum is the transverse size of the “patch” 
of a hadron over which the color is correlated. Thus, different emitter points separated by distances greater than Q~^ 
are uncorrelated in color. This random color orientation of different patches plays the role of the random phase of 
the individual sources in the previous discussion. 

In the next section we discuss how these features are reflected in the HBT correlations between emitted gluons. 


III. THE GLUON HBT 

As declared above, our goal is to understand gluon correlations. In this section we discuss the HBT signal for gluon 
emission in close analogy to the discussion of the previous section, highlighting its unique features. In the next section 
we will show that this discussion is closely paralleled by the glasma graph calculation, which will allow us to identify 
the HBT contribution, with all its distinguishing properties, in the current calculations. 


A. The gluon correlation function 

Like pions, gluons are emitted from the interaction region. However, gluon fields are real and consequently the 
emitter current is real in coordinate space (HBT of neutral pions would have to be treated similarly). Additionally, 
gluons carry an adjoint color index. Thus, in the same spirit as in the previous section, let us consider a superposition 
of N classical sources independently emitting gluons. Each source is translated to a different position Xi. The lack 
of coherence between the sources is encoded in a set of random adjoint unitary matrices {U\)ab, which determine the 
overall SU{Nc) “phase” of each source. 

The time and longitudinal coordinate dependence of all sources is identical, and resides in the factor (5(a;“'")5(x“). 
Therefore, hereafter the plus and minus light-cone components of all four-coordinates must bet set to 0. The exact 
delta function nature of the source is of course an approximation. In actual fact, even at very high energy the 
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gluon distribution in a hadronic wave function has a rapidity dependent structure. This is reflected in the Balitski- 
Fadin-Kuraev-Lipatov [17] (or Jalilian-Marian-Iancu-McLerran-Weigert-Leonidov-Kovner [18]) evolution of the gluon 
distributions with rapidity. However, the rate of the rapidity evolution is proportional to the strong coupling constant 
Us- Consequently, as long as the rapidity interval between the observed particles in the final state is smaller than l/os, 
approximate boost invariance holds and the source can be approximated by the delta function in the longitudinal 
direction. In the rest of this paper we will use this approximation. 

The emitter source function is 


N 


J,(p) = V (C/A),, Jo 6 (p) 


(3.1) 


with 


^ob(p) = '/ob(-p). (3.2) 

In this section, for simplicity of notation, we disregard the fact that the gluon field, and therefore also the gluon 
source, carry a transverse Lorentz (polarisation) index. This index will be restored in the next section. 

The momentum p in the expression Eq. (3.1) and thereafter is the two-component transverse momentum. The 
independence of the source function on the longitudinal momentum is an immediate consequence (upon Fourier 
transform) of the localisation of the source in time and longitudinal coordinate. For the same reason we have not 
allowed for any transverse motion of parts of the emitter. 

As before, we take (2.10) for the density matrix operator at a fixed value of the source J“ with the coherent state 
expression 


|J) =exp^i J cfpJaip) [ai(p)-f aa(-p)]^ | 0 ). (3.3) 

Here, the state |0) is the vacuum of the gluon Fock space. Since the source does not depend on the longitudinal 
momentum, the gluon creation and annihilation operators oA a here are integrated over rapidity [19]. In terms of the 
fundamental rapidity dependent gluon creation and annihilation operators ( 77 , p), one has (p) = J dr] da'^ ( 77 , p). 
This integration over rapidity is technically equivalent to approximating the source by a delta function in the longi¬ 
tudinal coordinate, and is valid with the same degree of accuracy. 

The immediate consequence of this is that both the single gluon spectrum and the inclusive two gluon spectrum 
do not depend on the rapidity of the emitted gluons. This is precisely the form of the density operator used in the 
glasma graph calculations of gluon correlations. 

The corresponding ensemble average of an arbitrary operator O is given as 

00 

tT{pd) ='^Pn 

iV=0 


Dn[J] {J[N; {x, Ua}]\0\J[N; {x, Ua}]). 


(3.4) 


For the case at hand, the measure is 


N 

Dn[J] = ( 3 - 5 ) 

with [dU] being the Haar measure, / [dU] = 1. 

One can calculate the single gluon correlator in a straightforward manner: 

CSO „ N 

(at(p)a 6 (k)) = ^ Piv / D^[J] (C/A)&d(C^l”')cae*'=-"-*P-"'Jo*,(p)Jod(k). (3.6) 

N—0 n,n' — l 

The integration over the unitary matrices can be performed by using the following formulae in the adjoint represen¬ 
tation: 


^ba^dc 


J [dU] {U)ab = 0, y [dU] {U)M{U^)ca 


m -I ■ 


(3.7) 
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Due to the orthogonality of the unitary matrices only the diagonal terms in the double sum over n and n' contribute. 
Thus, Eq.(3.6) reads 


oo « N 

(at(p)a,(k)) = Y.P^ 

= (3.8) 

where, as before, we have defined for convenience 

(J*(p)J(k)) = Ns j <fxnix)e-^^P-'^>^ J*M Joai^), (Ap)^'(k)) = Ng J d^xnix)e-^^P+'^>^ Joa(p) Joa(k).(3.9) 

For the single gluon spectrum, this gives 


(al(p)aa(p)) = (j;(p)Jc(p)). (3.10) 

An interesting and distinct property of the gluonic density operator, is that it yields a non-vanishing correlator of 
two-gluon annihilation operators, as well. This is a consequence of the reality of the gluon field (and source) in 
coordinate space. A simple calculation yields 

(aa(p)a&(k)) = - ^ (J(p)J(k)), (3.11) 

(at (p)at(k)) = - ^^^(A(p)A(k)). (3.12) 

The non-vanishing of these correlators is the consequence of the fact that the gluon operator does not carry any 
abelian conserved charge. Thus, the gluon density matrix, which in this case is a simple coherent state, does not 
possess a sharp ’’gluon number” but rather contains a superposition of Fock states with different numbers of gluons. 
For the two-particle spectrum we have 


CO p N 

(at(p)at(k)aa(p)a,(k)) = E E (UDaciUMUl^’)MUjr')d'b 

X e*P'(-n-.„de*fc(-™--^')j*^,(p)jg^(p)j*^,(k)j„,(k). (3.13) 

In the limit of very large number of incoherent sources the averaging over the SU{Nc) phases picks out three non¬ 
vanishing contributions, corresponding to {n = n'] m = m'), (n = m'\ m = n') and (n = to; n' = to'). The first two 
are similar to those discussed in the previous section, while the last one is new and corresponds to Wick contracting 
two creation operators and two annihilation operators: 


oo r { ^ 

(al(p)«l(k)aa(p)ab(k)) = Pn < W_d^Xi n{xi) 

AT _n I A—1 


1 


N 


N^O 
2 


{N^ - 1 )^ 


E {N^c - 1)^ ^Oa(P) W 


_|_ (]\r 2 _ 2 ^) ^i{p-k)-x„ ^-i{p-k)- 


n^m K 

Jo^(p)Joa(k) Jo%(k)Jo,(p) 


+ (iV2-l)e 


i{p+k)-Xn. g-iip+k)-x„ 


Jo\(p)Jo\(k) Job(k)Job(p) 


(3.14) 


Thus, the two-particle spectrum reads 

(iV(iV-l)) 


7^2 (P,k) = 


{Ny 


+ 


j(p) ) ( j(k) 

^ 'j*(p)j(k) 


m - 1 


A(p)A(k))(j(p)j(k) 


(3.15) 




























If the number of independent emitters is finite, an additional irreducible term is present in the spectrum: 


^2(p,k)= 


{N{N-1)) 

{Ny 


f / 


2 \ / 



^(p) 

)( 

J(k) 


irreducible 


^ Ns'^^ 

Defining the emission function 

S{x,K) = 


1 


m-l 


J*(p)J(k) 


J*(p)J*(k))(j(p)J(k) 


d'^y 


Ky 


2{2r. 


Ja{x + -y] Ja{x- -y 


we can write the normalized correlation function as 


C'(q,K) = l + 


1 


1 /(Px S{x,K) + / (PxS {x, I) e j (cc, e 


2iK-x 


m - 1 


J (PxS (x, K + ig) J (PxS (x, /’f — ig) 


(3.16) 


(3.17) 


(3.18) 


where we have neglected the irreducible term. As stated above, in all equations in this Subsection the plus and minus 
light-cone components of all coordinate vectors must bet set to zero. 

Eq. (3.18) differs from Eq. (2.28) in two aspects. Eirst, the number of gluons is — 1, and this is the origin of 
the suppression factor in front of the correlated term in Eq. (3.18). Second, gluon fields are real while the discussion 
in the previous section was for charged fields. As a result, Eq. (3.18) contains an extra correlation term - the second 
term in the numerator. As we will show below, this term contributes to the HBT enhancement between gluons with 
opposite transverse momenta and in this sense has a distinct signature. 


B. The transverse structure of the source and the peculiarities of the gluon HBT 

The transverse structure of the source Ja is what determines the HBT signal in Eq. (3.18). To understand it 
better, recall that the focus of our discussion is p-A collisions. The projectile proton carries into the collision region 
its distribution of classical gluon fields b\{x). This distribution in the proton wave function has some spatial size i?, 
which we identify with the gluonic size of the proton. Thus a typical gluon field configuration in the proton is a slowly 
varying function of the transverse coordinate, which varies on the scale R. 

The effect of the collision with the target nucleus is to rotate the field 6)j(x) {i is the transverse index) by an 
x-dependent SU{Nc) adjoint matrix U{x), so that immediately after the interaction the gluon field is Uab{x)h\{x). 
As we will see in the next section, in the CGC approach the role of the source J for the soft gluon emission is played 
to a good approximation by the gluon field at the instant following the interaction. Thus the source function that 
enters the calculation is closely related to the classical field, 

J*(x) « Uab{x) bl{x). (3.19) 

The averaging over the event ensemble in the calculation of the correlator (JJ) amounts to averaging over the 
distribution of the proton fields b{x) and also over the distribution of the unitary matrices U (x). The latter distribution 
is determined by the target wave function. Analogously to Eq. (3.9), we write 

{J{x)J{y)) = {{J{x)J{y))ti,))up) • (3.20) 

Assuming that the nucleus has a large saturation scale Qs ^ R~^ and that the eikonal scattering matrix U{x) 
typically varies on the spatial scale of we conclude that the source Ja is color correlated only on distance scales 
of order Qj^. The salient features of spatial structure of such source can be understood as visualizing the typical 
configuration Ja(x) as a collection of independent color sources, each having a fixed orientation in color space and 
the transverse size Qj^. The relative color orientation between the individual independent color sources in the event 
ensemble is completely random. The total area taken up by the source Ja(x) is R?^ thus the number of independent 
sources is Ns ^ Q^R^- 

Schematically (we will be more precise in the next section) we can encode the previous discussion in the following 
form of the source: 


N 

da{,x') — ^ ^ blab^l^doi^X ^i), 
i=l 


(3.21) 
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where ria is a fixed, unit length color vector common to all the individual sources, and the sources Jq{x) have a range 
of Q~^ <C R. The exact spatial dependence of the individual source Jo(x) is not important for our purposes. It affects 
the ’’irreducible” contribution to the correlation function, but has no effect on the HBT signal. As a simple model 
for the source distribution we can again use a Gaussian like in the previous section: 

_ 

n{xi) oc e (3.22) 


Note however that, as opposed to the previous section, the Gaussian profile in Eq. (3.22) involves only transverse 
coordinates. For the correlation function we then obtain 


C'(q,K) = l + 


1 






(3.23) 


The first term within the square bracket in this expression is the HBT signal of the “usual” kind. It leads to maximum 
of the correlation when the momenta of the two emitted gluons are almost equal to each other (with accuracy of R~^). 
The second term is somewhat unusual, since it maximises the correlation when the momenta of two emitted gluons 
are antiparallel to each other (with the same accuracy). This property of the gluon correlation function, namely that 
it is symmetric under the reversal of the direction of momentum of one of the gluons, has been discussed extensively 
in the literature, and is the consequence of the reality of the gluon production amplitude. This unusual term also 
constitutes part of the bona fide HBT correlation. 

To summarise this section we list the interesting properties of the gluon HBT correlation: 

• The correlation is long range in rapidity - it is equally strong when the rapidities of the two gluons are equal or 
when the difference between the two rapidities is large. 

• The correlation is symmetric under reversal of the direction of the transverse momentum of one of the gluons. 
Thus it is strongest when the transverse momenta of the two gluons are either parallel or antiparallel. 

• The HBT radius is of the order of the inverse gluonic size of the proton R~^. 

• The HBT signal dominates the correlation function (at small momentum difference) when the number of inco¬ 
herent emitters is large. Ns = Q^R^ ^ I- 


IV. THE HBT IN GLASMA GRAPHS 


In this section we identify the diagrams responsible for the HBT correlations in the ambient glasma graph calculation. 
Recall that in the eikonal approach the important contribution to the two gluon inclusive cross section in p-A collisions 
is given by 


{out\al^{p)al\k)al{p)al{k)\out) = 


X e 


dUW^[U] I DpW[p] (0|e*/^-'»(“)K^(“)+<(-")]at*(p)aJ^'(g)<(p)a^,(g) 

-»/x |Q^ 


(all coordinates and momenta in this Section are transverse and the corresponding integrals two-dimensional). Here, 
the state |0) is the vacuum of the soft-gluon Fock space in the presence of the ’’valence gluon” - a color charge 
density Pa{x), the Weizsacker-Williams field of the incoming projectile is given by bl^{x) = J d‘^y^{x,y)pa{y), and 
the color charge density configurations are distributed according to the probability density W\p\, while those of the 
color matrices in the target according to W'^\U]. The current J is given by 


Ja{x) = Uab{x)bl{x) - J d^y^{x, y)Uab{y)pb{y)- (4.2) 

Note that Eq. (4.1) is not the complete expression for inclusive two gluon production, but only the part where the 
two gluons are produced from two different Pomerons. This is the part of the production cross section which leads to 
the correlated result, and the only one that we will consider in this paper. 

The first observation is that the expression in Eq. (4.1) is exactly of the type discussed in the previous section. 
The averaging of the gluonic creation and annihilation operators is performed over a classical coherent state specified 
by the source Eq. (4.2), with subsequent averaging over the ensemble of sources distributed with some probability 
density W. The source is slightly different from the one we discussed in the previous section Eq. (3.21), but not 
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significantly so. The first term in Eq. (4.2) is precisely the source of Eq. (3.21). The second term is additional, 
however it is well known [20] that the emission of gluons with large transverse momentum (of the order of Qs) is 
dominated by the first term in Eq. (4.2). Thus, the presence of the additional term in Eq. (4.2) does not change 
significantly our expectations based on the discussion of the previous section. 

To calculate the two gluon inclusive cross section one has to average over the distribution of the color charge density, 
and also over the distribution of target helds Uab{x). The latter averaging introduces Q~^ as the typical transverse 
scale of variation of U(x). In current implementation one uses a Gaussian weight (the so called McLerran-Venugopalan 
model [21]) for averaging over the color charge density, and a similar factorizable distribution for the target average. 

After taking the expectation value over the soft gluon Hilbert space, the expression in Eq. (4.1) can be written as 

(outlal^(p)al^(k)a‘^(p)ai(k)lout} = (Jl‘(p)Ji^(k)J‘(p)Ji(k)}f^^[/} , (4.3) 

where the averaging over the distributions of p and U has to be performed. Using the factorizable model for averaging, 
this can be represented in terms of simple diagrams [7], which were used in [5] to calculate the gluon correlations 
numerically. The diagrams are depicted in Fig. 1. 



TYPE A1 



TYPE B 





TYPEC 


FIG. 1: Glasma graphs for the irreducible contribution (Types A1 and A2) and for HBT correlations (Types B and C). N{p — k) 
is proportional to the expectation value of the product of two U matrices, which has the meaning of the probability that the 
incoming gluon with transverse momentum k acquires transverse momentum p after scattering. 

It is easy to understand what is the role of the different diagrams in the framework that we have discussed in the 
last section (see also the discussions in [11]). First off, the diagrams that contain contributions of Eq. (3.18) are the 
ones shown in Fig. 1 named B and C. If in the calculation of the diagrams one assumes translational invariance of 
the projectile wave function, the loop integrals lead to momentum S functions: 

B oc (q — p), C oc (q + p). (4.4) 

Relaxing the translational invariance approximation, the delta functions are smeared over a scale of the size i? of 
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the projectile^. Thus diagram B represents the second term in Eq. (3.18), while diagram C the third term in Eq. 
(3.18). The identification of the diagrams B and C with the HBT correlation contribution is not new. It was known 
to the authors of [7] ([13]), and was also discussed subsequently in [14]. An earlier discussion of the diagrammatic 
representation of HBT can be found in [22]. The peculiarities of the gluon HBT signal was however not discussed in 
these references. 

It is also interesting to understand how the remaining diagrams, called Type A1 and A2, fit into the framework of 
our discussion. It was shown in [II] that the diagrams of Type A yield correlations due to Bose enhancement of gluons 
in the incoming wave function of the projectile and the target. This is the standard enhancement of the probability 
to find two incoming gluons in a hadronic wave function to have exactly the same transverse momentum. Due to 
momentum exchange during the scattering process, the enhancement is smeared over the momentum range of order 
Qs, but the correlation still remains. As discussed in [7] and [II], the contribution of these graphs to the correlation 
C(q,K) is suppressed by the factor l/{QlS), where S is the area of the smaller hadron (the proton). Identifying 
S with we see that this suppression factor is precisely the factor 1/Ns of the previous section. Thus the Type 
A contribution is identified with the irreducible term 'p^'reduciUe (3.16). Physically, this of course is clear, 

since the Type A contributions originate from correlations in the initial state, while the HBT correlations are due to 
independent uncorrelated emissions from the emitter formed immediately after the collision. 


V. DISCUSSION 

We have shown in this paper that the HBT signal in emission of gluons from the initial stages of p-A collisions has 
some interesting properties. 

It is boost invariant, that is, the HBT radius in the longitudinal direction is very small. In the boost invariant 
approximation used here, this radius is strictly zero and the correlation extends to arbitrary rapidity difference between 
the two emitted particles. Naturally, one expects that, at very large rapidity differences where the quantum evolution 
is important, the HBT signal must drop down. This should only happen for relatively large rapidity difference between 
the two gluons Nr] > 1/ag¬ 
in transverse space, the HBT radius is given by the gluonic radius of the proton. This, in principle, gives a possibility 
to directly measure the radius of the gluon distribution in the proton, assuming the signal is not washed out by final 
state effects. It is an interesting question what the gluonic radius of the proton is. The naive estimate would be of 
the order of the strong interaction radius ~ 0.8 — 0.9 fm. However there are many indications, particularly based on 
the DIS data, that the gluons in the proton are much more compact, taking up only about 0.3 fm (see e.g. [23]). A 
direct measurement of this quantity in p-A collisions would be extremely interesting. 

The HBT correlation is equally strong for two gluons with equal transverse momenta and momenta equal in 
magnitude but opposite in direction. This is a general feature of the gluon two-particle correlation. 

Finally, we have noted that the complete two gluon correlation within the glasma graph approach comprises two 
distinct contributions: the HBT and the irreducible initial state Bose enhancement. The relative importance of the 
two contributions is determined by the number of independent emitters Ns = . The two parts of the signal have 

different nature and different properties. The irreducible part is suppressed by 1/Ns, but leads to correlations whose 
width in the momentum space is determined by the saturation momentum Qg. On the other hand the HBT signal is 
unsuppressed, but is much narrower in momentum space, with the width R~^. The total signal that one expects is 
therefore a superposition of the two, and should have the general shape depicted on Fig. 2. It would be extremely 
interesting to see whether this double scale structure can be observed in the correlation. To do that, however, one 
would need to bin the signal into much narrower transverse momentum bin sizes than is currently done. 

We note that the separation of the two gluon correlation function into the HBT and the irreducible part only makes 
sense if the gluonic size of the proton is much larger than the inverse saturation momentum of the nucleus. Otherwise 


^ Note that the gluon fields in the proton can be correlated on a distance scale r < R. This corresponds to a ’’domain” picture where 
the proton contains several domains of color fields [8]. Nevertheless, as long as r 3> the scale that determines the HBT radius is 
the overall gluonic proton size R. Physically this is obvious since the only determining factor for HBT is the structure of the emitter, 
Eq. (3.21). In turn, the size of the correlated regions in the emitter is determined by the smallest spatial inhomogeneity scale which 
is . Mathematically, within the ’’glasma graph” calculation, the reason is that the momentum q — p in Eq. (4.4) is conjugate 

to the center of mass coordinate in the gluon correlator {a^{x)a{y)). The domain picture corresponds to the correlator of the type 
(at (x)a(y)} ~ fr{x — y)gR({x -|- y)/2), where the function fr vanishes when \x — y\ > r, while qr vanishes when \x -|- y\/2 > R. The 
smearing of the 5-function in the first relation in Eq. (4.4) clearly happens on the scale R. The same conclusion is reached for the 
second relation in Eq. (4.4). 
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FIG. 2: Schematic separation in q (or K) of the uncorrelated HBT (solid line) and irreducible initial-state Bose (dashed line) 
signals. Horizontal and vertical scales are arbitrary. 

Ns ~ 1 and the two signals become indistinguishable. Therefore one expects that this characterization is better suited 
for high multiplicity events, which correspond to the largest values of Qs- 

Of course, the usual word of caution is due here. All of the above discussion disregarded possible effects of final 
state interactions. Although the current thinking is predominantly that final state effects are very important in p-A 
collisions, nevertheless, as noted in the beginning, the ridge description based on initial state correlations is quite 
successful in reproducing available data [5]. In this vein one may hope that the final state effects do not completely 
obliterate the initial state signal, including the HBT correlations, at least for high enough transverse momentum. 
Admittedly, since the HBT signal is expected to be quite narrow, it is also very fragile and will be affected the most 
by any transverse momentum smearing in the final state. This includes final state interactions, but also fragmentation 
effects The fragility of the signal makes it rather tricky to predict in which kinematic regime it can be seen best. In 
particular, significantly increasing Qs also increases the number of produced particles and consequently the probability 
of final state interactions. On the other hand, for small Qs the difference between the HBT signal and the irreducible 
contribution to the correlation disappears. Thus, in order to see the signal it may be advantageous to study relatively 
forward production in p-A collisions in events with moderately high multiplicity. 
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of the saturation momentum, which made the HBT signal indistinguishable from the irreducible Bose enhancement correlation. The 
motivation for such smearing is precisely the effects of possible momentum transfer by final state interactions. 
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